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Abstract. The physical basis for interpreting observations of nebular morphology 
around massive stars in terms of the evolution of the central stars is reviewed, and 
examples are discussed, including NGC 6888, OMC-1, and r\ Carinae. 



1 Introduction 

The nebulae observed around massive, post-main sequence stars appear to con- 
sist of material ejected by the central stars during earlier phases of their evo- 
lution, rather than ambient interstellar matter. Models of these nebulae can be 
used to constrain the mass-loss history of the stars, giving an important input 
for stellar evolution models. Understanding the structure of these nebulae also 
clarifies the initial conditions for the resulting supernova remnants, which will 
interact with the circumstellar material for most of their observable lifetimes 
before encountering the surrounding interstellar medium. 

A strong stellar wind sweeps the surrounding interstellar gas into a stel- 
lar wind bubble as shown in Figure 1. The stellar wind expands freely until it 
reaches a termination shock. If this shock is adiabatic, the hot gas sweeps up 
the surrounding ISM into a dense shell, forming a pressure-driven or energy- 
conserving bubble that sweeps up the surrounding ISM into a dense shell grow- 
ing as R oc £ 3 / 5 in a uniform medium (Castor, McCray, & Weaver 1975). Should 
the termination shock be strongly radiative due to high densities or low ve- 
locities in the wind, the bubble only conserves momentum and will grow as 
R oc t 1 / 2 (Stei gman, Strittmatter, & Williams 1975). For more general discus- 
sions of blast waves in non-uniform media, see Ostriker & McKee (1988), and 
Bisnovatyi-Kogan & Silich (1995), as well as Koo & McKee (1990). 

When these stars leave the main sequence, they pass through phases of 
greatly increased mass loss. These slow, dense winds expand into the rarefied 
interior of the main sequence bubble until their ram pressure PwV^ drops below 
the pressure of the bubble. (Should the main sequence bubble have cooled rela- 
tively quickly, this may never occur.) As the mass loss rate and velocity of the 
central wind vary during the post-main sequence evolution of the central star, 
these denser winds can in turn be swept up, producing the observed ring nebulae 
around evolved massive stars. 

During their evolution, these nebulae are subject to a number of hydrody- 
namical instabilities, as well as thermal instabilities (e. g. Strickland & Blondin 
1995). I explain how understanding the physical basis of the hydrodynamical 
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Fig. 1. Stellar wind bubble structure. 



instabilities gives insight into the dynamics of observed nebulae. High-resolution 
observations of nebular morphology can thus be used to constrain the mass-loss 
history of the central star. 

I then describe how typical stellar evolutionary scenarioes can generate ob- 
served nebular morphologies, and show semi-analytic and numerical models de- 
rived from these scenarioes. For example, a star with a stellar wind varying from 
fast to slow and back again will have a clumpy circumstellar nebula due to hy- 
drodynamical instabilities in the shell (Garcfa-Segura, Langer & Mac Low 1996). 
Nonsphcrical winds and stellar motion can add to the morphological richness of 
the resulting nebulae, as in the nebula around rj Car (Langer, Garcia-Segura, & 
Mac Low 1998). A recent review of this topic is Frank (1998). 

2 Shell Instabilities 

Gas swept up by a stellar wind will usually be subject to different instabilities. 
An adiabatic, decelerating shell with a density contrast across the shock of less 
than 10 in a uniform medium is stable. However, relaxing any of these constraints 
will lead to instabilities as I now describe. 

2.1 Ray leigh- Taylor Instability 

If the swept-up shell is denser than the stellar wind, as will be true in virtually 
all cases where a shell exists at all, the shell will be subject to RT instabilities if 
the contact discontinuity between the shocked stellar wind and the shell accel- 
erates. This can be due to an external density gradient steeper than r~ 2 or to 
a sufficiently fast increase in the power of the central stellar wind, though these 
two mechanisms will lead to different shell morphologies, as I discuss below. 

The RT instability occurs when the effective gravity due to acceleration 
points from a denser to a more rarefied gas. We can understand its driving mech- 
anism by considering how the potential energy will change if we interchange a 
parcel of dense gas having mass mi with a parcel of more rarefied gas having 
mo < mi, as shown in Figure 2. The potential energy before the interchange 
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Fig. 2. Rayleigh- Taylor instability 

is given by = ui\gz + mog(z — Az) which is greater than the potential en- 
ergy after the interchange Ef — m n gz + m\g(z — Az) due to the difference in 
the masses. The decrease in potential energy drives an exponentially growing 
interchange of the two fluids. 

When a RT instability occurs due to an external density gradient, dense 
fragments of shell are left behind as the less dense interior expands out beyond 
them, creating the characteristic bubble and spike morphology seen, for example, 
in models of superbubble blowout (e. g. Mac Low, McCray, & Norman 1989). 
The Wolf-Rayet ring nebula NGC 6888 shown in Figure 3 provides another ex- 
ample. Here a fast, rarefied Wolf-Rayet wind has swept up the slow, dense red 
supergiant wind that preceded it. While it was still sweeping up the slow wind, 
it was marginally stable to RT instabilities. However, at the outer edge of the 
slow wind, the sharp density gradient triggers RT instabilities, as modelled by 
Garcfa-Segura & Mac Low (1995) with the astrophysical gas dynamics and mag- 
netohydrodynamics code ZEUS 1 (Stone & Norman 1992). 

On the other hand, when a RT instability occurs due to an increase in power 
of the driving wind, some of the dense fragments of shell actually get shot out 
ahead of the bulk of the fragmenting shell, producing a markedly different mor- 
phology (Stone, Xu, & Mundy 1995). Although these fragments represent only a 
small fraction of the total mass of the shell, they can produce a very striking set 
of bow shocks in their wake. An example of this occurring around one or more 
pre-main sequence stars is given by the "bullets" observed around OMC-1 (Lane 
1989, Allen & Burton 1993) in the Orion star forming region, as confirmed by 
McCaughrcan & Mac Low (1997). 

2.2 Vishniac Instability 

If a pressure-driven shell is decelerating, but thin, with a density contrast across 
the shock of at least 25 for a stellar-wind bubble expanding into a uniform 
medium (Ryu & Vishniac 1988), or 10 for a point explosion (Ryu & Vishniac 
1987), it will be subject to the Vishniac overstability (Vishniac 1983). This has 
been confirmed experimentally using blast waves generated by high-powered 
lasers propagating into gases with low and high adiabatic indexes (Grun et al. 
1991). 

1 Available by registration with the Laboratory for Computational Astrophysics at 
lca@ncsa.uiuc.edu 
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Fig. 3. Comparison between morphology of NGC 6888 in O [ill] (in an image taken by 
K. B. Kwitter with the Burrell-Schmidt telescope of the Warner and Swasey Observa- 
tory, Case Western Reserve University) , and a numerical simulation of RT instability 
due to a fast wind sweeping over the end of a slow, dense wind (Garcfa-Segura & Mac 
Low 1995). The model image shows a cross-section of the density structure in grayscale 
with black indicating high density and white low density. 

The mechanism of the Vishniac overstability can be understood by consid- 
ering a thin, decelerating shell driven from within by a high-pressure region, as 
shown in Figure 4. From within it is confined by thermal pressure acting normal 
to the shell surface, as adjacent regions can communicate with each other by 
sound waves, while from outside it is confined by ram pressure acting parallel 
to the velocity of propagation, as the shell moves supcrsonically into the sur- 
rounding gas. In equilibrium, these two forces remain in balance. Should the 
shell be perturbed, however, the thermal pressure will continue to act normally, 
but the ram pressure will now act obliquely, giving a transverse resultant force 
that drives material from "peaks" into "valleys" of the shell. The denser valleys 
will be decelerated less than the rarefied peaks, however, so that the positions 
of peaks and valleys are interchanged after some time. Vishniac (1983) showed 
that this overstable oscillation can grow as fast as i 1 / 2 . It saturates when the 
transverse flows in the shell become supersonic and form transverse shocks, so 
that the end result of Vishniac instability is a shell with transonic turbulence 
and moderate perturbations (Mac Low & Norman 1993). 

2.3 Nonlinear Thin Shell Instability 

Should the driving wind cool immediately behind its termination shock, for 
example because of exceptionally high mass-loss rates, it can form a deceler- 
ating shell that is momentum-driven rather than pressure-driven, so that it is 
effectively confined on both sides by ram pressure from shocks. Such a shell is 
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Fig. 4. Vishniac instability mechanism 

not subject to the Vishniac instability, and is, in fact, linearly stable. However 
Vishniac (1994) has shown that if the shell is strongly perturbed, it will still be 
subject to a nonlinear thin shell instability (NTSI). When the shocks are oblique 
enough to the direction of flow, they will bend the streamlines passing through 
them, so that mass is driven towards the extrema of the perturbation. Numerical 
simulations by Blondin & Marks (1996), using a piecewise parabolic hydrocode 
called VH-1, have shown that the end result is a catastrophic breakup of the 
shell into a turbulent layer that grows in time. 

3 A Final Example: Eta Carinae 

As an example of how knowledge of these different instabilities can be used to 
constrain the evolution of a star, consider the example of the Homunculus Neb- 
ula around r\ Car. Langer, Garcfa-Segura, & Mac Low (1998) computed several 
two-dimensional models of it using ZEUS, following a basic scenario in which 
a luminous blue star with a fast stellar wind undergoes an outburst during 
which it has a much slower and denser wind strongly shaped by rotation, as de- 
scribed by Bjorkman & Cassinelli (1993), but then reverts to its previous state 
with a fast, rarefied wind. They chose two different values for the post-outburst 
wind, one consistent with current observed values of M — 1.3 x 10~ 3 Mq yr _1 
and v w — 1300 km s _ , and one with a faster, lower mass loss wind having 
M = 1.7 x 1O~ 4 M yr" 1 and v w = 1800 km s -1 . As shown in Figure 5, the 
slower, denser wind cools upon shocking, forming a momentum-driven shell that 
fragments due to the NTSI, producing a sharp, spiky shell morphology. On the 
other hand, the faster wind does not cool completely, and forms a bubble subject 
to Vishniac instabilities, giving it a much more curved, cauliflower-like appear- 
ance. Comparison to the high-resolution observations (Humphreys & Davidson 
1994; Morse, Davidson, & Ebbets 1997) reveals that the actual morphology 
strongly resembles a three-dimensional version of the model with the faster wind. 
Langer et al. (1998) suggest that this reflects the typical behavior of the wind 
over the century since the outburst, and that the current wind properties are 
actually exceptional, and perhaps even indicative of another outburst on its way. 
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Fig. 5. Two-dimensional density distributions from models of r\ Car by Langer, 
Garcfa-Scgura, & Mac Low (1998) with faster and slower post-outburst wind show- 
ing Vishniac instabilities and the NTSI respectively. Note how the faster wind model 
resembles a cross-section through the cauliflower-like observed lobes. 



This suggestion is supported by the gradual brightening of r\ Car over the last 
decades (Humphreys & Davidson 1994). 

This work has made use of the NASA Astrophysical Data System Abstract 
Service. I thank the organizers for their invitation and their support of my at- 
tendence at this meeting. 
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